On the properties of steady states in turbulent axisymmetric flows 
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We experimentally study the properties of mean and most probable velocity fields in a turbulent 
von Karman flow. These fields are found to be described by two families of functions, as predicted 
by a recent statistical mechanics study of 3D axisymmetric flows. We show that these functions 
depend on the viscosity and on the forcing. Furthermore, when the Reynolds number is increased, 
we exhibit a tendency for Beltramization of the flow, i.e. a velocity-vorticity alignment. This result 
provides a first experimental evidence of nonlinearity depletion in non- homogeneous non-isotropic 
turbulent flow. 
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Introduction A yet unanswered question in statistical 
physics is whether stationary out-of-equilibrium systems 
share any resemblance with classical equilibrium systems. 
A good paradigm to explore this question is offered by 
turbulent flows. Incompressible flows subject to statis- 
tically stationary forcing generally reach a steady state, 
in a statistical sense, independent of the initial condi- 
tions. Most of the current turbulence modeling is de- 
voted to the understanding of this state, and the role 
of the velocity fluctuations or velocity smallest scales in 
its construction. A distinguished feature of stationary 
large Reynolds number turbulent flows is the presence of 
coherent structures, under the shape of vortices in 2D 
P, 0| , or vorticity thin tubes in 3D j^] . These structures 
correspond to regions where vorticity is locally almost 
aligned with velocity — phenomenon referred to as Bel- 
tramization 0, 0. This tendency to velocity-vorticity 
alignment induces depletion of nonlinearities in Navier- 
Stokes equations. Interestingly enough, similar depletion 
of nonlinearity is also observed in the inviscid counter- 
part of the Navier-Stokes equation — the so-called Euler 
equation — resulting in a slowing down of the vorticity 
blow-up with respect to rigorous estimates 0. A theo- 
retical question of interest is therefore whether turbulent 
flows have a natural tendency for nonlinearity depiction, 
and whether this is a characteristic of the steady states. 

The answer to the first issue is ambiguous. On one 
hand, a tendency for Beltramization has indeed been 
observed in some numerical simulations of stationary, 
nearly isotropic turbulence 0, On the other hand, 
similar study performed on more general flows — such as 
boundary layer Q — provided little evidence of such a 
property .. A similar conclusion has been reached using 
experimental data |lfl | , with slightly less reliability owing 
to the difficulty to measure vorticity. As for the second 
issue, it has only been partially explored in special ge- 
ometries. In 2D, equilibrium states of the Euler equations 
have been classified through statistical mechanics princi- 
ple by Robert and his collaborators [ill ll2T ] . They indeed 
correspond to Beltrami solution. In 3D axisymmetric 
flows — an intermediate situation between 2D and 3D — a 



similar task has been undertaken by Leprovost et al. |13| . 
In the ideal case — force free, inviscid — they proved the 
existence of an infinite number of conserved quantities, in 
addition to the energy and the helicity. They also showed 
the existence of an infinite number of equilibrium states, 
depending on the conserved quantities. Among them, 
a Beltrami state is obtained, when only the energy and 
the helicity are conserved. They postulated that even in 
the presence of small but finite viscosity and forcing , 
this feature remains valid, with selection of steady states 
through boundary conditions and forcing. 

The purpose of the present Letter is to perform an ex- 
perimental check of these findings. We find that mean 
and most probable velocity fields of axisymmetric flows 
can indeed be characterized by two families of functions, 
as predicted in Ref. 0. These functions depend on the 
viscosity and forcing. As the Reynolds number is in- 
creased, they evolve towards the functions corresponding 
to a Beltrami state, providing a first evidence of non- 
linearity depletion in a non-homogeneous, non-isotropic 
system. 

Theoretical background and definitions Consider an 
incompressible axisymmetric flow, with velocity compo- 
nents in a cylindrical referential (v r , vg, v z ). Due to in- 
compressibility and axisymmetry, only two functions are 
sufficient to describe the flow, namely the angular mo- 
mentum — o~(r, z) = rvg — and either, the stream func- 
tion ^(r,z), such that (v r ,v z ) = — V x ^eg, or the az- 
imuthal vorticity LUg(r, z). Using a variational method, 
Leprovost et al. showed that the steady states solution 
of axisymmetric Euler equations, in the force-free case 
obey: 

FF' 

o- = F{V); £- — = G(V); with £ = u e /r, (1) 

where F and G are arbitrary functions linked with con- 
servation laws of the system — Casimirs of <r, generalized 
helicity 01 . For example, the steady state correspond- 
ing to conservation of energy and helicity is such that F 
is linear and G = 0, resulting in a linear relation between 
vorticity and velocity (v = Xlo). This is a Beltrami flow, 
given by : 
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where A and m are free parameters and Jo and J\ are 
first order Bessel functions. 

Considering further the thermodynamics of the sys- 
tem, Leprovost et al. proved that the equilibrium states 
at some fixed coarse-grained scale are such that the most 
probable velocity flow {mpvf) is a stationary state of the 
Euler equation, i.e. satisfies Eq.(QJ. Finally, they pos- 
tulated that these force-free, inviscid results can be ex- 
tended to the stationary states in the presence of both 
viscosity and forcing, providing the function F and G are 
selected through boundary conditions and forcing. 

Experimental setup In order to check the theoretical 
predictions, we have worked with a simple "axisymmet- 
ric" configuration: the von Karman flow generated by two 
counter-rotating impellers in a cylindrical vessel. The 
cylinder radius and height are respectively R = 100 mm 
and H = 180 mm (distance between inner faces of im- 
pellers). The impellers consist of 185 mm diameter disks 
fitted with sixteen 20 mm high curved blades. More 
details about the experimental setup can be found in 
Rcf. 0. The impellers rotation frequencies are both 
set equal to / to get exact counter-rotating regime. We 
define two forcing associated with the concave (resp. con- 
vex) face of the blades going forward, denoted in the se- 
quel by direction (— ) (resp. (+) ). The working fluid is 
either water or glycerol at different dilution rates. The 
resulting accessible Reynolds numbers (Re = 2irfR 2 v~ 1 
with v the kinematic viscosity) vary from 10 2 to 3 x 10 5 . 
In the exact counter-rotating regime, whatever the im- 
pellers, the flow is divided into two toric cells separated 
by an azimuthal shear layer. This setup is invariant un- 
der rotations of n (TZ^) around any radial axis passing 
through the center of the cylinder. The time-averaged ve- 
locity fields we consider hereafter are T^-invariant and 
axisymmetric |l5j . Velocity measurements are done with 
a DANTEC Laser Doppler Velocimetry (LDV) system. 

Data processing The LDV data only provide the ax- 
ial and azimuthal velocity components on a 170 points 
grid covering half a meridian plane, through time series 
of about 200,000 randomly sampled values at each grid 
point. From this time series, it is straightforward to get 
time-averaged axial and azimuthal velocities. The re- 
maining radial component of the mean velocity field can 
then be obtained using the incompressibility and axisym- 
metry: this procedure has been later validated through 
direct measurements of the radial velocity with a Parti- 
cles Image Velocimetry (PIV) system on the same flow. 
We also use the time series to extract at each point his- 
tograms of the axial and azimuthal velocity, and compute 



the flow of most probable velocities (fmpv) |r| > obtained 
by taking, at each point, the most probable value of the 
velocity. The v r component is once more derived from the 
two other components using continuity equation. Fig- 
ure n compares the two flows computed from the same 
data set: they have the same overall structure, but dif- 
fer through the size of the middle shear layer, which is 
thinner in the case of the fmpv. This is due to the fact 
that the chaotic wandering of the shear layer around its 
average position are significantly less probable as one con- 
siders positions further away from the equatorial plane. 
From a theoretical point of view, the two flows differ in 
the sense that the mean velocity is not a solution of the 
Navier-Stokes or Euler equations — because of the fluctu- 
ations, which induce a Reynolds stress, especially in the 
shear layer — , while the mpvf is a stationary solution of 
the Euler equations. Therefore, one can expect theoriti- 
cal predictions regarding the structure of the stationary 
state to be less accurate in the case of the mean flow, as 
the Reynolds number and the fluctuations increase. 
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FIG. 1: Velocity profile for Re = 2.5 x 10 5 , direction (-): col- 
ored contour for azimuthal component, arrow representation 
for poloidal part. Top: fmpv. Bottom: time-averaged flow. 

We first compute the F function by looking at a as a 
function of 'J, and then, we use this estimate to evaluate 
G, by looking at £2 = £ — as a function of This 
procedure is likely to induce a lot of noise in the estimate 
of G. To check its robustness, we have tested it on a 
Beltrami flow (Eq. [5J with a superimposed level of noise 
comparable to the level of fluctuations in the flow (see 
Figure |2J). Even in the presence of noise the fits give the 
correct shapes: F — A '5 and G = 0. 

After this test, we apply this procedure on the real 
data, with a further correction, motivated by the follow- 
ing remark. Figure displays two typical plot of F : one 
over the whole apparatus, and one over a 50% portion of 
the flow obtained after removing regions close to the im- 
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FIG. 2: a and £2 = £ — ^S- versus * for the Beltrami flow 
defined by equation with a superimposed white noise of 
amplitude 60%. We have represented the F (resp. G) fit on 
the left (resp. right) figure. 

pellers and along the vessel. These parts correspond to 
locations where the viscosity and the forcing (neglected 
in Ref. are principally at work, and larger devia- 

tions from theoretical predictions can be expected. One 
sees that while data on the whole vessel display signifi- 
cant scatter, preventing the outcome of a well-defined F, 
the restricted data gather onto a cubic-shaped function 
fitted by a two parameters cubic: F(^) = pi^ + p3 1 F 3 . 
This fit is then used to obtain G. In the sequel, we obtain 
F and G through fits over this 50% portion of the flow 
away from the boundaries. 



similarly. As Reynolds number increases, the F cubic 
curves collapse on the Beltrami line. 

At Re — 2.5 x 10 5 , it is not even possible to distin- 
guish the F fit from a straight line and this for the two 
directions of rotation. In the case of G, there is a differ- 
ence between the two cases. For the fmpv, the best fit 
is linear, with a slope smaller than the noise level, this 
is consistent with G = 0. For the time-averaged case, 
the shape of G is very different. All the fits present a 
plateau around zero. The greater the Reynolds number, 
the wider this plateau, so the wider the range of W - 
or equivalently the volume of considered flow — where G 
is very close to zero, i.e., very close to the Beltrami G 
function. 

The difference in the shape of G can be explained as 
follows : according to equation QJ, G(W) is the difference 
of two terms. One linked to F is necessarily a cubic, 
while the second, connected to the azimuthal vorticity is 
almost a noise in this area of the flow. In the fmpv case, 
the magnitude of the cubic term is small compared to 
the vorticity, and we only see the noise coming from the 
vorticity. In the time-averaged case, the cubic is much 
larger, so that it determines the behavior of G. 
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FIG. 3: a vs "if for experimental time-averaged flow in direc- 
tion (+) at Re = 2100. Left: for the whole flow. Right: for 
r < 0.81, —0.56 < z < 0.56, corresponding to 50% of the flow 
volume. The remaining points are clearly gathering along a 
cubic-shaped function F. 

Results We present results computed for both time- 
averaged flow and fmpv for the two directions of rotation, 
at different Reynolds numbers. Figure 0] presents the F 
fits obtained in each case. We have included on each 
graph a straight line tangent to experimental curves for 
\f r = as a reference to a Beltrami flow. The experi- 
mental data are much more scattered in the fmpv case, 
debasing the fit accuracy compared to the mean flow — 
R 2 ~ 0.94 instead of 0.99. 

In figure we present the G fits obtained in the same 
conditions. The best fit of the fmpv G functions is lin- 
ear whereas cubic-shaped functions are required for time 
averaged G. The plots of G consist of wide noisy bands 
surrounding the fits in both cases. The noise for G is 
larger than for F, since it is a result of a two-step proce- 
dure including a spatial derivation. 

Considering the variation with Reynolds number of F, 
both the fmpv case and the time-averaged case behave 




FIG. 4: F fit. (a) and (b): fmpv direction ( — ) and (+) respec- 
tively, (c) and (d): time-averaged field direction ( — ) and (+) 
respectively. Legend for all figures: Re = 100 (thin ••■), Re = 

150 (thin ■ —), Re = 2100 (thin ), Re = 9100 (thin —), 

Re = 2.5 x 10 5 (thick — ). The thin dotted straight line is 
a Beltrami. As F is odd, we display it for positive values of 
$ only. Note that the fit at Re = 9100 for fmpv case in di- 
rection (+) could not be obtained from the data due to large 
scattering. 

A second observation is the dependence of F and G 
on the forcing. Comparison of Fig. |U-(a,c) and|3J(b,d), 
obtained for two different forcing, shows that the slope 
of F slightly depends on the forcing (pi — 5 for (— ) 
direction, p\ — 4 for (+) direction). For any forcing, G 
remains close to zero at the measurement scale. To check 
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FIG. 5: G fit. (a) and (b): fmpv direction ( — ) and (+) re- 
spectively, c and d: time-averaged field direction (— ) and (+) 
respectively. Same legend as Fig. ^] 

this dependence further, we have conducted an additional 
experiment at high Reynolds number, with much smaller 
impellers. As shown in Figure we obtained a different 
shape for F, with a different inflection. 

Discussion and perspectives An inviscid force-free 
theory developed by Leprovost et al. predicted a charac- 
terization of the steady state in axisymmetric turbulent 
flows through two functions F and G. fn our experi- 
ments, we have confirmed this, and measured these func- 
tions for different forcing, and different viscosities, using 
two different fields as diagnostic: the fmpv as predicted 
and more surprisingly the time-averaged field, even in a 
region where Reynolds stresses are not negligible. This 
could be due to the quasi Gaussian distributions of ve- 
locities. 

First, we have observed that the functions are well- 
defined only in the portion of the flow remote from the 
boundaries, where forcing and dissipation take place. 
The steady states of the von Karman flow can be de- 
scribed by a suitable statistical analysis of equivalent 
equilibrium states (e.g. assuming zero viscosity and no 
forcing, as in Ref. |ig). Nevertheless, we have seen 
that forcing and dissipation do influence the steady state 
regime through the selection of the characteristic func- 
tions. Specifically, we have observed that the forcing se- 
lects the specific shape of the functions (linear, cubic, ...), 
while dissipation acts as a sort of self-similar zoom and se- 
lect the portion of the curve actually explored by the flow. 
This suggests that stationary out-of-equilibrium systems 
like turbulent flows are universal in a weaker sense than 
in ordinary equilibrium systems : they can be described 
in a universal manner through general functions (like F 
and G) determining the "equation of state". However, 
these functions are non-universal since they depend on 
the fine details of the system (dissipation and forcing). 



Finally, we have shown that the evolution for increas- 
ing Reynolds number is towards a Beltrami state, with 
depletion of nonlinearities. This evolution is more ob- 
vious when considering the fmpv, as expected from the 
thermodynamics analysis developed by Leprovost et al.. 
To our knowledge, this is the first experimental evidence 
of nonlincarity depletion in a non-homogeneous, non- 
isotropic turbulence, where it is very challenging to mea- 
sure simultaneously all components of velocity and vor- 
ticity. Our use of the functions characterizing the steady 
states enables us to lower the experimental constraints 
for such a check. 
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FIG. 6: F fit for 100mm diameter impellers fitted with 
straight blades. Left: fmpv, Right: time-averaged. Re = 
1.9 x 10 5 (■ ■ ■), Re = 2.6 x 10 5 ( — ), Re = 5 x 10 5 (— ). 
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